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Abstract. We discuss the effects of strong magnetic fields through Landau quantization of 
electrons on the structure and stability of nuclei in neutron star crust. In strong magnetic 
fields, this leads to the enhancement of the electron number density with respect to the zero 
field case. We obtain the sequence of equilibrium nuclei of the outer crust in the presence of 
strong magnetic fields adopting most recent versions of the experimental and theoretical nuclear 
mass tables. For B ~ 10^^ G, it is found that some new nuclei appear in the sequence and some 
nuclei disappear from the sequence compared with the zero field case. 

Further we investigate the stability of nuclei in the inner crust in the presence of strong 
magnetic fields using the Thomas- Fermi model. The coexistence of two phases of nuclear matter 
- liquid and gas, is considered in this case. The proton number density is significantly enhanced 
in strong magnetic fields B ^ lO^^G through the charge neutrality. We find nuclei with larger 
mass number in the presence of strong magnetic fields than those of the zero field. These results 
might have important implications for the transport properties of the crust in magnetars. 



1. Introduction 

The discovery of a new class of neutron stars with very strong magnetic fields called magnetars 
has greatly enhanced the interest in the study of neutron star properties in the presence of strong 
magnetic fields [T]. Their surface magnetic fields could be > 10^^, as predicted by observations 
on soft gamma-ray repeaters and anomalous x-ray pulsars [21 [3] . Such strong magnetic fields 
might be generated by dynamo processes in newly born neutron star In the core region the 
magnetic field may be even higher, the limiting value (Bmax) is obtained by the scalar virial 
theorem |5J. For a typical neutron star (M = 1.5Mq, R = 15 km) this value is B^ax ~ 10^^ 
G. Such high magnetic fields can have significant effects on the equilibrium nuclear composition 
and equation of state in neutron star crust and interior [5l [6] . 

Nonmagnetic equilibrium composition and equation of state for the outer crust was reported 
in a seminal paper by Baym, Pethick and Sutherland (BPS) Outer crust contains nuclei 
arranged in a body-centered (bcc) lattice immersed in a gas of free electrons which are relativistic 
above the density p ~ 10 ''g cm~^. Though the lattice effect is small on the equation of state, 
it changes the equilibrium nucleus to a heavier one and lowers the total energy of the system 
by reducing the coulomb energy of the nucleus. At p ~ 10^ g cm~^, ^^Fe is present as the 
equilibrium nucleus, but with increasing density equilibrium nuclei become more and more 
neutron rich through electron capture process. At a density p ~ 4 x 10^^ g cm~^ neutrons begin 
to drip out of nuclei - this is called neutron drip point. The inner crust begins from here. In the 
inner crust nuclei are immersed in a neutron gas as well as a uniform background of electrons. 



Nuclei are also arranged in a lattice in the inner crust. The composition and the equation of 
state of the inner crust were earlier calculated by different groups [U [9] . 

The composition and equation of state of the outer crust of nonaccreting cold neutron stars in 
the presence of strong magnetic fields were first studied by Lai and Shapiro [5j . In the presence 
of a magnetic field the motion of electrons perpendicular to the field get quantized into Landau 
orbitals. This causes the electron density to change which in turn modifies the coulomb energy. 
If the magnetic field is very strong then electrons occupy only the low-lying Landau levels and 
it may affect the sequence of nuclei and the equation of state as well as any nonequilibrium 
/3-processes [S]. However, there is no calculation of the inner crust composition and equation of 
state in the presence of magnetic fields. 

This paper is organised in the following way. We revisit the magnetic BPS [S] adopting 
recent experimental and theoretical nuclear mass tables in Section II. The inner crust in strong 
magnetic fields is discussed in Section III. We conclude in Section IV. 



2. Magnetic BPS Model 

We revisit the BPS model to find the sequence of equilibrium nuclei and calculate the equation 
state of the outer crust in the presence of strong magnetic fields B ~ 10^^ G [5j. In this 
calculation, we include the finite size effect in the lattice energy and adopt recent experimental 
and theoretical mass tables. Nuclei are arranged in a bcc lattice in the outer crust. The Wigner- 
Seitz (WS) approximation is adopted in this calculation. Each lattice volume is replaced by a 
spherical cell and contains one nucleus at the center. Each cell is taken to be charge neutral 
such that Z number of electrons presents in it, where Z is the nuclear charge. The Coulomb 
interaction between cells is neglected. To find an equilibrium nucleus (A,Z) at a given pressure 
P one has to minimize the Gibbs free energy per nucleon with respect to A and Z. The total 
energy density is given by 

Etot = riNiWN + Wl) + Seine) . (1) 
The energy of the nucleus (including rest mass energy of nucleons) is 

Wn = mn{A - Z) + rUpZ - bA , (2) 

where h is the binding energy per nucleon. Experimental nuclear masses are obtained from the 
atomic mass table compiled by Audi, Wapstra and Thibault (2003) [10]. For the rest of nuclei 
we use the theoretical extrapolation of Moller et al (1995) [Hj. Wl is the lattice energy of the 
cell and is given by 

M.,.-jL^ri4f^rv (3) 
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Here rc is the cell radius and ~ r^A^/^ (ro ~1.16 fm) is the nuclear radius. The first term in 
Wl is the lattice energy for point nuclei and the second term is the correction due to the finite 
size of the nucleus (assuming a uniform proton charge distribution in the nucleus). Further Eg 
is the electron energy density and P is the total pressure given by 

P = Pe + ^WlUn- (4) 

The nucleon number density tt-at is related to the baryon number density rif, as 

Ub = AuN, (5) 

and the charge neutrality condition gives 



Ue = ZriM- (6) 




Figure 1. Electron number density is plotted with mass density for different magnetic field 
strengths. 



In the presence of a magnetic field, the electron motion is Landau quantized in the plane 
perpendicular to the field. We take the magnetic field (^) along Z-direction and assume 
that it is uniform throughout the outer crust. If the field strength exceeds a critical value 
Be = m1/e ~ 4.414 X lO^^G, then electrons become relativistic. The energy eigenvalue of 
relativistic electrons in a quantizing magnetic field is given by 

Ee{i^,Pz) = [pI + mlil + 2uB,)\ , (7) 

where Pz is the Z-component of momentum, = B/Bc, and v is the Landau quantum number. 
The number density of electrons in a magnetic field is calculated as 



max 



where the spin degeneracy g^^ = 1 for = and = 2 for all other Landau levels. The 
maximum Landau quantum number Vmax is given by 

where jie is the electron chemical potential. Energy density of electrons is given by, 

^<^ = ^ ^ 9u Ee{v,Pz)dpz , (10) 

and the pressure of the electron gas is 

Pe = flene - £e ■ (H) 

At a fixed pressure P, we minimize the Gibbs free energy per nucleon 

_ Etot + P _ Wn + 4./3WL + Zfle 



varying A and Z of a nucleus. 

Now we present our results for several values of magnetic fields = B / Be = 0, 100, 1000 
and 10000 where B^ = 4.414 x lO^^G. In Fig. 1, electron number density is plotted with mass 
density for the above mentioned values of magnetic fields. For B^, < 1000, large numbers of 
Landau levels are populated. Consequently there is no significant change in the electron number 
density compared with that of the field free case. But there is a jump in the electron number 
density for B^: = case when only zeroth Landau level is populated at ~ 2 x 10^^ g/cm^. On 
the other hand, for B^ = 10^, only zeroth Landau level is populated over the whole mass density 
range. In this case, the phase space modification due to the Landau quantization leads to the 
strong enhancement of electron number density with respect to the zero field case. 





Figure 2. Proton number is shown as a function of neutron number for different theoretical 
nuclear models with and without magnetic fields. 

We obtain the sequence of equilibrium nuclei minimizing the Gibbs free energy per nucleon 
for S = and B = lO^^G. These results are obtained with and without the lattice energy 
correction. We compare our findings at B^ = 10^ with those of zero field case and find that 
some new equilibrium nuclei such as IgSr and Ig^Pd appear and some nuclei such as ^^Ni and 
^^Ni disappear in the presence of the magnetic field. Further we note that the maximum density 
{Pmax) upto which an equilibrium nucleus can exist, increases as the field strength increases. 
As a result the neutron drip point is shifted from 4.34 x 10^^ g/cm^ in zero field to 4.92 x 10^^ 
g/cm^ in the strong magnetic field. The lattice energy correction influences our results in strong 
magnetic fields. It is worth mentioning here that our results are different from those of earlier 
calculation [5] because we have adopted most recent experimental and theoretical nuclear mass 
tables. For S = 0, we note that ijg^Sr which was present in the previous calculation [5] is absent 
in our case. A comparison of our results at zero field with those of Lai and Shapiro [5] shows that 
equilibrium nuclei have higher Pmax in the former case. Further we performed our calculation 
at higher magnetic fields than the previous calculation [5]. 

Figure 2 displays the proton number as a function of neutron number. We compare sequences 
of equilibrium nuclei calculated with experimental nuclear masses from atomic mass table of 
Audi, Wapstra and Thibault [10] and theoretical nuclear mass models for those nuclei which 
are not listed in the experimental mass table. We use the theoretical model of Moller et al. 
|11] and relativistic mean field model with NL3 set [121 [T3] . It is evident from the figure that 



for zero magnetic field our calculation of equilibrium nuclei initially agrees with those of the 
relativistic model calculation because nuclear masses are obtained from the experimental mass 
table. However, both calculations for zero magnetic field differ considerably beyond N=50 due 
to differences in theoretical mass tables used here. 
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Figure 3. Mass number is displayed as a function of average baryon density. 



3. Inner Crust in Magnetic Fields 

Next we investigate the stability of nuclei in the inner crust in the presence of strong magnetic 
fields using the Thomas-Fermi model. In this case nuclei are immersed in a neutron gas as well 
as in a uniform background of electrons. The coexistence of liquid and gas phases of nuclear 
matter is considered here. In the inner crust, nuclei are also arranged in a lattice. The Wigner- 
Seitz approximation is adopted in this calculation. However, the spherical cell in which protons 
and neutrons coexist does not define a nucleus. The nucleus is realised after subtraction of the 
gas part from the cell. This is done following the procedure of Bonche, Levit and Vautherin [14j. 
In this case, the density profile of liquid plus gas system and that of the gas are obtained in a 
self-consistent way and the liquid part which ultimately becomes the nucleus, is calculated as 
the difference of two solutions. We assume /^-equilibrium = A^p + t^e- Electrons are affected 
by quantizing magnetic fields like the outer crust. The proton density in the cell is affected by 
magnetic fields through the charge neutrality condition. We minimise the free energy which is 
a function of average baryon density (n^) and proton fraction (Yp), in the cell 

J'inb, Yp) = J[n + ec + Eeldr , (13) 

where H is nuclear energy density functional, Ec is the coulomb energy density and Eg is the 
energy density of electrons. The nuclear energy density is calculated using SKM* nucleon-nucleon 
interaction [15]. We perform this calculation for temperature T=0. 

Now we present our results related to the inner crust. We fix the average baryon density rib 
and achieve the /3-equilibrium changing the proton fraction Yp. Next we find the minimum of 
the free energy varying the cell size. We perform this calculation for B = and = 10^, 10^. 
We do not find any appreciable difference in our results for B = and -B* = 10'^. However, the 
cell size is appreciably reduced in case of B^ = 10^ . In this case only zeroth Landau level is 




populated below the density 0.06 fm~^. This leads to the large enhancement of proton fraction 
in the cell through the Landau quantization of electrons. Above this density, electrons populate 
finite number of Landau levels and the cell size and proton fraction in the magnetic field case 
become equal to those of the zero field. The size of the cell as well as proton fraction decrease 
as density increases. 

Figure 3 displays mass number (Aceii) of nuclear cluster as a function of average baryon 
density for S = and = 10^, 10^. These are obtained after minimising the free energy at 
each density point. It is noted that at any density point the mass number of the nuclear cluster 
in magnetic field = 10^ is lower than the corresponding nuclear cluster at zero magnetic field. 
Further we obtain the liquid part of the cell after subtraction of the gas part and define this as 
a nucleus. We find nuclei with larger mass number in the presence of a strong magnetic field 
~ 10^''G compared with the zero field case. Further nuclei are more bound in strong magnetic 
fields than the corresponding nuclei in i? = at the same baryon density. 

4. Summary and Conclusions 

We have revisited the BPS model of outer crust in the presence of strong magnetic fields ~ 10^® G 
or more using the recent experimental mass table. Further we have included the correction in 
the lattice energy due to the finite size of a nucleus. For zero magnetic field case, it is noted that 
maximum densities for heavier nuclei in this calculation are higher than those of the previous 
calculation [5]. In the presence of a strong magnetic field, there are modifications in the sequence 
of nuclei compared with the zero field case. We have further investigated the inner crust in the 
presence of strong mag netic fields. For B = lO^^G, it has been observed that proton fraction 
is enhanced at lower densities and mass numbers of nuclear clusters after subtracting the gas 
part are higher than those of the zero field. It would be worth investigating the implications of 
our results for the transport properties such as thermal and electrical conductivities and shear 
viscosity of the crust in magnetars. 
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